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Abstract 
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peculiar values of M and N with new tricritical points. When adding weak disorder, the 
results we obtain tend to show that disorder makes the models decouple. Therefore, no 
relations emerges, at a perturbation level, between for example the disordered q\ x (72- 
state Potts model and the two disordered q\, (72-state Potts models (q± / (72), despite 
their central charges are similar according to recent numerical investigations. 
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1 Introduction 



In order to understand the role of impurities and inhomogeneities in real physics systems, 
many statistical models with quenched randomness have been proposed. Therefore, the 
effect of randomness of continuous phase transitions have been of great interest for many 
years. The main question is to determine if the randomness leaves unchanged the critical 
properties of the pure system. One prediction concerning models with random bonds comes 
from the Harris conjecture [I]], which states that bond randomness changes the values of 
the critical exponents only if the specific-heat exponent a is positive. The effects of bond 
randomness has been studied intensively in the 2D Ising model (The Harris criterion does 
not enable to be conclusive about the relevance of disorder) first by Dotsenko and Dotsenko 
||. They represent the near critical point by a Gross-Neveu model and found an asymp- 
totic behavior for the specific heat in Zn(ln(4r) with \t\ the reduced temperature. This has 
been confirmed numerically 0. For 2D random Potts model (the Harris criterion predicts 
here that randomness is relevant and changes the critical behavior), conformal field theory 
techniques have been a powerful tool to compute in pertubation the shifts in the critical 
exponents j|, |5|. The numerical results have proved this approach to be accurate (see [Q 
and for more recent results). All these results concern simple models having a second 
order phase transition in the pure case. 

However, it has been realized recently, that a weak bond randomness can also have strong 
effects on 2D systems possessing a first order transition. Hence, following earlier work of 
Imry and Wortis ||, Hui and Berker have shown with a phenomenological renormalisation 
group that bond randomness will induce a second order transition (a vanishing of the latent 
heat) in a system that would have undergone a first order one |§. This result has been stated 
on a more rigorous level by Aizenman and Wehr [10|. At this level, a question naturally 
arises concerning the universality class of this second order transition. In order to test these 
predictions, Cardy has studied a system presenting a fluctuation driven first order transition, 
namely a transition which is expected to be continuous at the mean field level, but become 
first order when fluctuations are included. Therefore, Cardy was able to show that the ad- 
dition of weak randomnes on a system of iV coupled Ising models makes the system flow to 
N decoupled Ising models [|ITJ. This result has been extended to the case of iV (N > 2) 
coupled Potts models by pujol |12| (the case N = 2 is quite peculiar and deserves a special 
attention fl2| , p!3|). The result was a non-Ising like second order transition. The universality 
class depends on the sign of the coupling between the models and therefore these systems 
seem to violate the Imry- Wortis argument ||. 

Naturally, the q— state Potts model with q > 4 appears more suitable to test and study 
the effects of disorder, because the first order transition is by now of mean field type. The 
hint is that there is now analytical approach able to control the effect of disorder except at 
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large q where a mapping to the random field Ising model shows the absence of any latent heat 
in 2D 0. For q not too large, only numerical simulations are able to check the Imry-Wortis 
arguments. Hence, Chen et al. jnj] have investigated the 8— state Potts model using Monte 
Carlo simulations and confirmed the transition to be continuous but also found numerical 
values of the critical exponents consistent with those of the Ising model. Nevertheless, recent 
numerical studies of Cardy and Jacobsen and of Picco |15| are in a clear disagreement 
with the latter conclusion. The magnetization exponent ^ is found to vary continuously 
with q. Moreover, it has also been shown, using the powerful tool of conformal invariance 
combined with finite size scaling, that the values of the central charges of disordered q— state 
Potts models are related to one another by a factorization law c(q = q\ x q 2 ) = c(q\) + c(q 2 ) 
15, [7j. Therefore, the numerical measure of central charges can not distinguish between 



a disordered q = q\ x g 2 -state Potts model and two decoupled q\ and q 2 disordered Potts 
models. As the only analytical results concern the behavior of similar coupled q— state Potts 
model (q = 2, 3) under weak disorder, it should be interesting to extend these results to 
different disordered Potts models in order to compare the critical behavior of coupled q\ and 
q^-state Potts models under weak randomness with the disordered qi x g 2 -state Potts model. 
In order to bring some answers to this issue, the main goal of this paper is therefore the study 
of the critical behavior of M Ising models coupled to N q— state Potts models (g = 3,4) 
with disorder. Nevertheless, a special attention will be first paid to the pure case. 
Indeed, it seems important to understand in a general context the critical behavior of different 
minimal models notably coupled by their energy density. This question was recently adressed 
by Leclair et al. |16[ (see also [17]]) in the context of integrable perturbations theories. Namely 
could the mixing of two or more critical models create a new critical behavior ? Moreover, 
it might be a way to approach the critical behavior of some 3D systems. For example, when 
we couple one Ising model to iV g-state Potts model, we obtain using one loop calculations 



a non-trivial flow with new tricritical points |L8|. Nevertheless, the most puzzling situation 
concerns the unusual critical behavior of one Ising model coupled to an XY model, currently 



referred as Ising-XY in the literature [19]. This model is expected to include in its range of 
parameters a whole class of systems with Z 2 xU(l) symmetry like the 2D fully frustrated XY 
model [^U| the 2D arrays of josephson junctions the ANNNXY model and so on. The 
Ising-XY model have in its phase diagram a continuous critical line with simultaneous XY 
and Ising ordering where the critical exponent associated to the Ising magnetization varies 
continuously and differs significantly from the pure Ising one. Therefore, along this critical 
line, we have two models at their critical temperature coupled by their energy density and the 
numerical results indicate an unusual behavior. Those examples shows that an energy-energy 
coupling term may have drastic effects in certain situations. Therefore, before introducing 
weak disorder, we carefully study the critical behavior of M Ising models to iV q— state Potts 
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models. We thus confirm, using two loops calculations, the existence of the two tricritical 
points quoted above. These results also extend when we consider M 3— state Potts models 
coupled to iV 4— state Potts models. 

The paper is written in a self-contained way and completes the work presented in fT3| , pl| . It is 
organized as follows. In section 2, we analyse the critical behavior of M Ising models coupled 
to N Potts models. A specific attention is paid to the cases M = 1,2. Near their critical 
points, these models are presented by perturbed conformal theory. Then, in section 3, we 
study the effects of weak randomness on these models. We use the replica trick method. We 
show that disorder makes the models decouple. Therefore, our analysis confirms no apparent 
relations (at least perturbatively!) between the disordered q% x g 2 -state Potts model and the 
corresponding coupled qi, q 2 disordered Potts model, despite their central charge are similar. 
Section 4 contains a summary of the results and some future perspectives. All the integrals 
involved in the computation of the beta functions, the renormalisations of the spin and 
energy operators are worked out explicitly in the appendices. 

2 Coupled minimal models without disorder 
2.1 The model 

The model we indend to study in this section consists of a superposition of M Ising and 
N g-state Potts models coupled by their energy density. The system possess a Zif x 
symmetry. The Hamiltonian has the following form 




where c, d belong to [1, • ■ ■ , M + N]. H° and H\ represent respectively one pure Ising model 
and one pure q— state Potts model at their respective critical temperature. e c corresponds 
to the energy operator of the c th model (either Ising or Potts depending on the value of c). 
Therefore, the partition function reads as 

m m+n . 

m+n -J2 S o,a~ H So,b+ / Y,9cd,e c {z)e d {z)d 2 z . 

Z =U Tr a , s , a II T n, Si>h e a=1 b = M+1 ^ (2.2) 

o=l b=M+l 

This model can be described as a pure conformal field theory perturbed by the interaction 
term which is quadratic in the e operators. The infrared behavior of the coupling constants 
of the complete model (g^ and also the different masses that have not been included in Q2.1|) ) 
can be analysed perturbatively. In a similar way, a correlation function like < O(0)O(R) >, 
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where O is some local operator, can be expanded perturbatively like: 

< O(0)O{R) >=< O(0)O(R) > + < SfO^OiR) > +^ < SjO(0)O(R) > + • • • 
where <>o means the expectation value taken with respect to the bare action and 

Si = j Hi{z)d 2 z = J Y,9cdE c (z)e\z)d 2 z , (2.3) 

c^=d 

the interaction term. 

2.2 Coulomb gas formulation 

In order to compute the integrals involved in the calculus of correlations functions, we used 
the Coulomb-gas representation of a conformal field theory ||24j| . In this representation, the 
central charge is written as c = \ + e, where e will be used as a short distance regulator 
for the integrals involved in correlation functions calculations. In our calculations, we are 
obliged to introduce two UV regulators ei,e2, since we mix two different minimal models. 
In addition, we also use an I.R. cut -off r which is useful to derive the RG equations. Then 
the limit e x — ► and a finite value for e 2 corresponds to Ising models coupled to q— state 
Potts models. The integrals will thus be expressed as finite series in ei,e2 with coefficients 
depending on r. 

In the Coulomb gas representation, the central charge c will be characterized in the following 
by the parameter o? + = ^rj — f + e with 

c = 1 - 24«q ; a± = a ± \J a§ + 1; (2.4) 
= —1. 

For the pure Ising model a 2 + = | and c = \ while for the 3— states Potts model a\ — |, 
c = | and e = — t=. The 4— states Potts model can be considered as the limiting case of 

5 15 ° 

this perturbative scheme ||, it has c=l and e = — |. The vertex operators are defined by 
V nm {x) = e 10 ™" 1 ®^ where $(x) is a free scalar field and the a nm are defined by 

OL nm = ^(1 - n)a^ + ^(1 - m)a + . (2.5) 
The conformal dimension of the operator V nm (x) is A nm = —otnfnOt nm with 

anrn = 2tt ~ "nm = ^(1 + Ti)OL- + ^(1 + m)a + . (2.6) 

The spin field a can be represented by the vertex operator V P}P -i whereas corresponds to 
the energy operator e. Note that in the Ising case the a operator could also be represented 
by the V^i- Therefore, the spin operator can be represented in the general case (Ising or 
Potts) by V k ,k-i where k = ^±^. We thus have A = 2 for V21 and A = \ for V p>p -\. All 
these notations generalise in a straightforward way when we mix two minimal models since 
we introduce two UV regulators E\ and e 2 . 
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2.3 Renormalisation group equations 

In this section, we deal with the renormalisation of the <jj and Si (i = 1,2) operators. 
Nevertheless, we first need to compute the renormalisation of the couplings constants gij. 

2.3.1 Beta functions 

When we have only one coupling constant go, its renormalisation is determined directly by 
a perturbative computation, g is also given by the O.A. producing 

go J e a (z)e b (z)d 2 z + ^ (g J e a (z)e b (z)d 2 z^j H = g J e a (z)e b (z)d 2 z , (2.7) 

with g = go + A 2 g 2 + • • ■ where A 2 comes from the contraction 



i J e a (z)e b (z)d 2 z J e c (z)e d (z)d 2 z = A 2 J e a (z)e b (z)d 2 



(2.8) 



This procedure generalizes to two loops and more. In this article, we restrict to two loops 
calculations. 

All this formulation can be extended in a straighforward way to our case. We suppose that 
9cd = 9i for c,d G [1,M], g cd = g 2 for c,d e [M + 1,M + N] and finally g cd = g 12 for 
c G [l..M],d G [M + 1,M + N] or vice versa. The technical details are explained in the 
appendix A. We find 



gi(r) = r 



-3ei I „0 



g° 1 -MM-2)(g° 1 



,0\2 



-3ei 



+ 8tt 2 (M-2)(^) 3 ' 



36! 
-6ei 



-^N(g° 12 



\2' 



-3e 2 



3e- 



3ei 



+ 32vr 2 iV(M-2)^(^ 2 ) 2 



1 + 



2{M-2) 



3ei 



^-3(ei+£2) r I 



3(ei + e 2 ; 



3ei + 3e 2 



(2.9) 



+ 167T 2 A^(^ 2 



-3(ei+e 2 ) 

3(ei + e 2 ) 



2 

3ei 



+ 16vr 2 iV(iV - l)^ 



,0/^.0 \2' 



-6e2 



12; 9e 2 



g 2 (r) = r-^(g 2 -An(N-2)(g° 2 ) 



r ~3e2 r — 3ei 



0\3 



+ 8vr 2 (iV-2)(^) 



r 



3e 2 



1 + 



2{N-2) 



3e 2 

r -3(ei+£2) 

+ 32 7r 2 M(iV-2)g 2 (g° 2 ) 2 

<j(ei + e 2 j 

„-3(ei+e2) r O 



1 1 

3ei 3e 2 



+ lQrr 2 M(M - l)^) 2 - 



(2.10) 



-6ei 



9e 2 
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-3(ei+e 2 ) [ _ 



-3ei 



-3e 2 



g» 2 -4n(M-l)g»g» 2 - 



-4^-1)^ 



+ 4n 2 (M-l)(g°)y i2 



+ ^(N-l)(g°yg 12 



3ei 
3e 2 



1 + 
1 + 



3ei 
2(M- 1 



3e 2 



36! 

2(iV- 1) 



+ 167r 2 (iV-l)(M-l)^ 2 V 2 



3e 2 

„— 3(ei+£2) 

3(ei + e 2 ) 



1 



1 



3ei 3e 2 



(2.11) 



lA^0N3 ^ 3(£1+£2) r 1 I 
12 > 3( Cl + 6 2 



+ lQir 2 (N -l)(M -l)(g 

r -3(ei+£2) 



3ei 3e 2 



+ 8vr 2 (,? 2 ) 3 



3 ei + e 2 ) 



», *r „ 2(M-1) 2(JV-i; 

M + iV-2 + — - + 



3ei 3e 2 

We have multiplied each equations by a power of r in order to have dimensionless coupling 
constants ft(r). From these expressions, it is then easy to obtain the three beta fonctions 
associated to g± : g 2 , gvi- We obtain the following result 

^91 , _ , fit n \ J i Ar J C\/r n\ „3 at „ „2 , i „3\ 



ft 
&2 



dilnl) 



e 1 gi + (M — 2) gf + N g{ 2 - (M - 2) g{ - N g ig f 2 + o(^), (2.12) 
e 2 92 + {N- 2) g\ + 4ttM £ 2 2 - (N - 2) (? 3 - M ^ 2 2 + o(^ 3 ), (2.13) 



rf(/n/) = ^ 12 + (M ~ 1} ^ 12 + (iV " 1} ^ 12 " \ {M + N ~ 2) ^ 2 



~(M - 1) ^ 12 - i(iV - 1) glg l2 + oQ? 3 ) . 



(2.14) 



In the above equations, we have made for simplicity the changements g — > Aug and e, — > — 3ej. 
If we consider the case of M Ising models coupled to iV 3— state Potts models, we take the 
limits e 1 — > and e 2 — > f . The fixed point structure is studied in a further section. 



2.3.2 Renormalisation of 

Under the action of the perturbation term, the operators will be renormalised into e^. It 
means that we have to compute the 2x2 matrix Z £ defined by e' = Z £ e, where e' are two 
components vectors. Indeed, since the interaction term mixes different energy operators, the 
matrix Z £ contains off-diagonal terms. In order to compute these matrix function, we add 
source terms to e in the action, namely mass terms. Therefore, the renormalisation of the 
operator e is equivalent to the renormalisation of the mass term fa defined by 

fa E = m° Z £ E, (2.15) 

with 



(. M r M+N \ 

/E^)A/ E e b 2 (z)d 2 z) 
J o=l b=M+l J 
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As for the computation of the beta fonctions, we can compute in perturbation the renor- 
malisation of m. The details of the computation are presented in appendix B. We thus 
have 



-3ei 



TOl (r)r- 1+ i £l = m Al-kK{M-l)gf- + 4vr 2 (M - 1)(^ ) 2 ^— 

\ 3ei v 3ei 



-wi r AM -6 
1 + 



3d 



+ 8vr 2 iV(^ 2 



\2 



r -3(ei+£2) 

3(e! + e 2 ) 



2M 2M-1 
1 + + 

3ei 3e 2 



(2.16) 



/ r ~3e2 r -3(ei+€2) 

+«2 -4^2-— + 8tt 2 N(M - l)g° ig 12 — 

V 3e 2 9exe 2 



-6e2 ' 



+ 8ir 2 N(N - 



9e 



wi2 r r 



-l+|e2 



-3e 2 



mo 



47T(iV-l)^ - 



3e 2 



+ 4vr 2 (iV- l)(g 



0\2 ' 



-6e2 



3e s 



1 + 



4iV-6 



3e 2 



\2 
12; 



-3(ei+e 2 ) 



3(ei + e 2 ) 



+ra? —AuMg 



-3ei 



12" 



3e a 



2iV 2 7V- 
1 + + — 

3e 2 3e x 



8n 2 M(N - l)g»g» 2 



(2.17) 



^-3(ei+e2) 

9eie 2 



8vr 2 M(M - 1)^ 2 



9e 2 



Here again, we have multiplied the two renormalised masses by a power of r in order to 
obtain a dimensionless coupling constant. From these equations and the renormalisations of 
the coupling constants (|2.9| , [2.10| , [2.11|) , we obtain the RG equations at two loops associated 
with the matrix Z e (we have again used the normalisations g — > 47rg) 



d\og{Z £ ) / ( M -l) gi -\{M-l)gl-\Ng 



12 



Ng 



12 



log(r) 



Mg 



12 



(N-l) gi -\{M-l)gl-\Mg 



12 



(2.18) 

The logarithm of the matrix has been defined by its series expansion. We clearly see that 
both energy operators of the two models are mixed under the perturbation term. Therefore, 
the bare operators £j are not the true energy eigen-operators of the problem. 



2.3.3 Renormalisation of a { 

In a similar way, the spin operators Oi will be renormalised in a\. We can thus define the 
fonctions Z ai , Z a ^ by a\ = Z ai Ui with i = 1,2 (the interaction term does not mix the spin 
operators). A convenient way to compute these functions is to introduce two magnetic field 
h® which couple to the spins Then, we have to calculate the effect of the perturbation 
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M+N 

term J2 9cd I d 2 z e c (z)e d (z) on the coupling terms J2 J2 <Ji i {z)d 2 z. More precisely, the 

c,d i=l,2 ai=l 

corrections to the spin operators will come from the operator algebra between the e and 
a operators. Therefore, as in the previous subsections, we can compute these corrections 
perturbatively in power of gij up to third order. The computation is detailed in the appendix 
C. We find 



^(rJr-T-^O = h\ l + (M-l)(#°) 2 7r 2 ' 



1 + _(2_A 1 )^ 



e 2 + e 1 (l + ~(2-A 1 )ft) 



-12(M-l)(M-2)(^ 



-8iV(iV - 1)7^)^) 



\N{M - iy{g Q 12 ) 2 gl 



r 



-9e x 



9ex 

^-3ei-6e2 

3ei + 6e 2 

„— 6ei— 3e2 



6ei + 3e 2 

and in a similar way by exchanging (1, M) <-> (2, N) 



l + -(2-Ai)W 
"l + |-(l + ^(2-A 1 )^ 

:i + -)(i + |(2-a 1 )^) 

e 2 3 



:(1 



(2.19) 



62- 



/i 2 (r)r-ir- a ^ = /i° 1 + (N - 1)G? 2 )V 



3M, 2 2 r ~ 3 ( £1+£2 ) 



12 J 



71 



3(e! + e 2 : 



l + -(2-A 2 )ft 



ei + e a (l + -(2- \ 2 )K) 



-12(N - 1)(N - 2)(g° 2 ) 3 7T 3 



,-9e2 r 



-8M(M - 1)tt^ 2 )^) 



r 



9e 2 

5ei— 3e2 



l + -(2-A 2 )W 



M(N-l)7t 3 (g° 12 ) 2 g° 2 



6ei + 3e 2 

— 3ei — 6e2 r 



3ei + 6e 2 



l + ^(l + -(2-A 2 )^) 

;i + -)(i + |(2-a 2 )^)- 

ei 3 



(2.20) 



where we have defined 



ft 



r2 H) r2 H)' 

As usual, the multiplicative term r~^~ a ^^ in front of hi(r) is introduced in order to deal 
with dimensionless parameters. The function a(ej) is a function of depending on which 
representation of the spin field we take in the Coulomb gas picture (see section 2.2). Now, 
we deduce from the above equations and from ( |2.9| . p.l0| . p.ll|) , the RG equations associated 



to Z ai , Z a2 (with ei 
dlog(Z ai 



d\og(r) 



-3e f ) 

: (M-l)^(r)7r 2 ei 



i + -(2-x 1 )n 
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+Ng 



TV 



12V 



e 2 + ei (l + -(2-A 1 )^) 



+4(M - 1)(M - 2)<K 6 gl{r) + 4iV(iV - l)<K 6 g{ 2 {r)g 2 {r) (2.21) 
+4iY(M-l)7r 3 ^ 2 (r)^(r) 



dlog(Z, 



0"2 / 



dlog(r) 



(N-l)g 2 2 (r)n 2 e 2 



l + -(2-\ 2 )K 



+Mg 2 12 (r 



71- 



e 1 + e 2 (l + -(2-\ 2 )K) 



+4(iV - 1)(N - 2)7T 3 gf(r) + AM(M - 1)tt 3 g{ 2 {r)gi{r 



3^2 



(2.22) 



-4M(iV - l)w 3 g 2 12 (r)g 2 (r 



2.4 Fixed point structure 

In this section, we analyse the fixed points structure associated with the beta fonctions 
defined in fl2.12 - 2T4|) . We focus in this paper on the case of Ising models coupled to 3 or 



4-states Potts models. Indeed, the 3— state Potts model corresponds to e = | in our reduced 
normalisation whereas the 4— state Potts model corresponds to e = 1. One can wonder about 
the validity of the perturbation expansion for the latter case. In fact, the beta fonctions have 
still fixed points order by order. And surprisingly, the approximation still remains quite good 
qualitatively. For example, if we compare the magnetic exponent of the disordered 4-state 
Potts model predicted by the theory with numerical results, we find an error less than 5% 



We now take the limit e\ — > and e 2 — >• e in the equations ( [2.12 - 2~T4] ). They reduce to 
A = ^ = (M -2) g 2 1 + N g\ 2 ~ (M -2) g\- N g l9 2 2 , 

Pi 



d{lnl) 
dg X2 



eg 2 + (N-2) g 2 + 4ttM g{ 2 - (N - 2) g s 2 - M g 2 g 2 12 , 



(2.23) 



^ = lfrn = ^gi2 + {M-l)g igi2 + {N-l)g 2 g l2 -\{M + N-2) g\ 2 
a{lnl) 2 2 

-~(M - 1) g\g X2 -~(N-1) g 2 2 g 12 . 
We consider in the following different peculiar values of (M,N). 



2.4.1 Case M = 1 

It corresponds to one Ising model coupled to one or several Potts model. 

• Let us first analyse the case of one Ising model couple to one g-state Potts model. Then, 
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the equation ( |2.14j ) reduces to 



dg 



12 



12 



d(lnl) 



2 9i2 + o{ 



A) 



(2.24) 



We clearly see that the system is then driven in a strong coupling regime indicating probably 
a first order transition though a non-perturbative fixed point cannot be ruled out. A similar 



situation has been encountered when one couple two g-state Potts model fl2|| . Moreover, in 
the latter case, it can be proved exactly that we have a mass gap generation |T7j. Conse- 
quently, when one superposes a critical Ising model to a critical 3 or 4 states potts model, no 
new critical behaviour is found in pertubation when coupling both models by their energy 
density. Of course, new non-perturbative fixed points cannot be excluded by this method. 
Therefore, It suggests (perturbatively) that no line with simultaneous disordering of the Ising 
and Potts order parameter survives in the phase diagram of this model contrary to those of 
the FFXY model |]I9|. Similar conclusions can be drawn with two different Potts models. 



• We now consider the case of one Ising model coupled to N > 1 Potts models. We take 
M = 1 in the equations (|2.23|) . We find the following fixed points: 



9i= 92- 

* * 

9l = 9u 
9*i = 



- 9n = 
= 0; g* 2 = 



+ 



2(7V-1) 



(N-2) (iV-2) 2 

Tie 2 (3iV-l) 
8(N-l) 2 



if n > 2 



(2.25) 
(2.26) 



a* - -£ 4- ( N + 1 ^ 2 
^2 2(N-1) ~ 8(iV-l) 2 



(2.27) 



9l2 



r 2 e 2 (3jV-l) 
8%/iV(jV-l) 2 



2{N-1) 

with Ti, r 2 = ±1. There are four new fixed points (|2.27|) . We can study the stability of each of 
these fixed points by re-expressing the solutions of ( |2.23| ) around these fixed points. Namely, 
we write gi = g* + 5gi and keep the smallest order in e. This gives a linear system 5g = A 5g, 
with 5g a three component vector and A a real 3x3 matrix. Therefore, the informations 
about the stability of each fixed points is obtained by calculating the eigenvalues of A in each 
cases. By the way, we can see that the trivial fixed point (|2.25|) is unstable, and that ( |2.26| ) 
if stable in the (g 2 , gn) plane and marginal in the gi direction. It corresponds to the infrared 
fixed point of N coupled Potts models. Concerning the other four fixed points ( p.27| ), two 
are unstable in two directions (when t% = —1) and the order two (ri = +1) are hyperbolic 
fixed points, namely unstable only in one direction. They constitute new tricritical points. 
The local stable plane is defined by the vectors 



r 2 \ 2n(n -1+1-1 

1,0,0) ; (0, V ' ' ,1) 

y/n(n - 1) 
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in the (gi,g 2 ,gi 2 ) space. A projection of the flow in the (gi 2 , #2) plane has been given in 
Figure 1 for the case N = 3. The flow is symmetric along the #2 axis. If initially g® > then 
the flow is always driven in a strong coupling regime indicating mass generation as for the 
simple iV-colored Potts model Nevertheless, if g® < and small, the trajectory will be 
first attracted by one of the tri critical points 7\ or T 2 (depending on the sign of g® 2 ) an d then 
flow toward the stable fixed point of the iV-colored Potts model. The case iV = 2 is quite 
different since the fixed point ( |2.26| ) is rejected to the infinity, therefore the system always 
flow in a strong coupling regime except on the stable surface in the (gi, g 2 , #12 ) space. It 
is a reminiscence of the exactly integrable 2-colored Potts model which is always massive [17|. 



If we consider the case of one 3-state Potts model coupled to N 4-state models, we can 
derive exactly similar conclusions. The four tricritical points are now defined by (at first 
order in ei, e 2 ) 



91 



TlJ4 + 



N 



(ei + e 2 )(JV(e 2 -ei) + 2ei) 



9l ~ 2(N-1) 



(2.28) 



9h 



2(N-1) 



e 1 + e 2 ))(iV(e 2 -e 1 )+2e 1 ) 



We find a similar structure of fixed points as the previous case, namely two are stable in one 
direction, two are stable in two directions. 

In order to better characterise these fixed points, we can compute the corrections to the 
critical exponents of the pure case. The energy operators are completly mixed at the tricrit- 
ical points and are therefore not the most appropriate variables, especially for a numerical 
simulation. Nevertheless, we can compute the renormalisation of the critical exponents of 
the spin-spin correlation functions. Hence, we use the Callan-Symanzik equations which give 
the form of the correlation functions. For one coupling constant g, we have 



la a (9) 



2 J /3(a) 

< ai(0)ai(sL) > a , go = e 9 ° s~ 2 ^ < a % (0)(j % (L) > r>9(s) 

where we used the notation : dl ^f^ = Ja{g) ', g{a) = go and g(s) is defined by 



dg 



(2.29) 



dln(r) 



/ P(g)dg = log(s) ; r = sa, 



with a the lattice cut-off scale. L is an arbitrary scale which can be fixed to one lattice 
spacing a of a true statistical model. The dependence in s of the term < <7i(0)<7j(L) > r , g (s) 
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is thus negligible and can be considered as a constant. We want to compute the corrections 

to the critical exponents close to the infrared tricritical points defined in (|2.27|) . Hence, the 

g(s) , ) 

integral is dominated by the region g ~ g* and we have / ^jW-dg ~ 7 S (g*) log(s). We thus 

go {a> 

obtain 

< a(0)a(s) > ff0 ~ s -( 2A - 2 7.(9*)) (2.30) 

Therefore, — 2 r y CTi (g*) corresponds to the correction of the magnetic critical exponent of the 
operator <Tj. Using the equations ( |2.21| , j2.22| ) with e\ — 0, e 2 = e, M — 1, we find 



7«iGf) = iVye(^ 2 ) 2 + 4iV(iV-l)7r 3 (^ 2 ) 2 ^ = + O(e 4 ) (2.31) 

Tff2 (<f) = (iV-l)(^)V e [l + 2^] + (^ 2 ) 2 y6[l + 4^] 

+4(iV - 1)(AT - 2)tt 3 (^) 3 + 4(iV - 1)tt 3 (^ 2 ) 2 ^ 

~ 128(iV-l) 2 + 32{N-lf'^ + U{e h {I - 6I) 

with e = | for the 3-state Potts model and e = 1 for the 4-state Potts model. We thus find 
the surprising result that there is no renormalisation of the critical exponent A (T1 associated 
to the Ising spin at third order in e! Yet, the corrections to A a2 are non zero and different 



from those of the infrared fixed point of iV coupled Potts models fl2|. Therefore, it sems 
that the physics at these tricritical points is non-trivial. It should be interesting to analyse 
and test numerically these predictions despite the corrections to the magnetization are very 
small. Before considering the other cases, we mention one remark. Suppose we have a the- 
ory described by one of the tricritical fixed points Tj, i — 1,2 (see Figure 1), then a small 
perturbation can send the flow to the fixed point D characterized by decoupled Potts and 
Ising models. The central charge at this point Co is thus the sum of the central charges of 



each model. Therefore according to Zamolodchikov c-theorem p3| , the central charge at Tj 
(ctJ verifies > cp. A similar scenario could be imagined for the FFXY model with a new 
tricritical point inducing a strong cross-over regime. Such a scenario would justify strong 



finite size effects and cffxy > f as it has been established numerically [19|, |20| . 



2.4.2 Case M > 1 

We analyse in this subsection different peculiar cases and the general one. 
• Let first begin with M = 2, N = 1. The RG equations reduce in this case to 



n dgi 2 2 
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^ 2 = djtl) = ^2-^ + 2^ + ^-2^ (2.33) 

dg 12 e 1 3 1 2 

^ = d{M) = 2 9l2 + 9l9l2 -2^-2 9 ' 912 - 



The fixed points are 



e + e 2 



g m i =?;<£ = < Q ; = o (2.34) 



We recover the line of fixed point parametrised by g\ which corresponds to the Ashkin- Teller 



model model |TT|]. We can study (perturbatively) the stability of this line of fixed point 
especially according to the variable (around g\ = e + e 2 corresponding to the fixed point 
of the Potts model). We have 

^fH| + Sl *-i( S ;)^ P.35) 



We thus see that the line of fixed point is stable only for g\ < — § + We have drawn 
on Figure 2 a shematic representation of the flow in the (^1,^12) plane). The point P(<7* = 
— I + y\ 92 = e + e2 5 9*2 = 0) separates the line parametrized by g\ in two parts, a stable one 
and an unstable one. At the fixed point P, we can compute the corrections to the spin and 
energy critical exponents. We find on one hand no correction to the Ising magnetization and 
on the other hand 

2A'(P) = 2A £ - 7ei (P) = 2A £1 + e + o(e 2 ). 



This result may be checked numerically. 

• If we now consider the case M = 2, N = 2, we can draw similar conclusions than 
the previous case except that the fixed point g\ = e + e 2 has disappeared and the line 
parametrized by g* is therefore unstable in the gi direction. 

• The case M = 2, N > 2 can be studied in a similar way. We find the following fixed 
point 

gl =? ; gl = I ^ ; 9u = o (2.36) 

2 

The value g^ = -^F^ + (N-2) 2 corres P on ds to the fixed point of the iV colored Potts model 
0. We can show that the latter is stable for g% < otherwise we have a strong coupling 
regime indicating probably mass generation. We recover the Askhin- Teller fixed line which 
is stable for gi < independently of gu. 
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• Finally, if we consider the case M > 2. We find the fixed points g\ = g\ 2 = 0. The 
value of g 2 depends on the value of N we consider. 

When we consider only one Potts model, we find that g\ 2 = is unstable whereas g\ = is 
stable for g® < 0. 

For N = 2, the only fixed point is the trivial and unstable one. Therefore, the flow is always 
driven in a strong coupling regime. 

For N > 2, it depends on the initial conditions we consider. For g® < and g 2 < the 

2 

fixed points g\ = 0; g\ = + nv-2) 2 > #12 = is stable. Nevertheless, the situation is unin- 
teresting because the Ising and Potts models are completly decorrelated in the infrared limit. 



We have thus studied in this general section the critical behavior of M Ising models and 
iV Potts model coupled each other by their energy density. New non-trivial situations have 
been exhibited for M = 1 and M — 2. We refer to the Figures 1 and 2. The case of one 
Ising model coupled to several potts models is interesting in so far as it could offer a possible 
scenario for the puzzling physics of the Ising- XY model |B5 . 



We have also seen that in many cases, the flow is driven in a strong coupling regime (it 
often depends on initial conditions). It should be interesting to analyse the effects of weak 
disorder on such systems in order to test the Imry-Wortis arguments. It will be the matter 
of the next section. 



3 Coupled minimal models with disorder 

In this section, we analyse the effects of weak disorder on the previous systems, namely M 
gi-state Potts models coupled to N g 2 -state Potts models, with qi,q2 E {2,3,4}. We will 
pay special attention to the case q% = 2; q 2 = 3. 



3.1 The model 

The Hamiltonian of the pure system have the following form (see also (|2.1|) ) 




M . M+N . 

+m x 53 / d 2 xel + m 2 ^ d 2 xe b 2 (3.1) 

0=1 J b=M+l J 

The masses rrii correspond to the reduced temperatures. The addition of a position depen- 
dent random coupling constant is equivalent to consider position dependent random mass 
terms m; — > rrii(x) with rrii(x) = mi and rrii(x)mj(y) = Aij5(x — y). A^- represents the 



15 



2x2 symmetric covariance matrix whose elements are strictly positive as it should be. If 
we have considered a diagonal covariance matrix corresponding to independent disorders for 
each models, then Ai 2 would have been generated by the R.G. equations. We then apply the 
replicated method by introducing n copies of the system and averaging Gaussian distribu- 
tions for irti{x). We finally obtain the sum of nM q± and nN g 2 -state Potts models coupled 
by their energy densities. The Hamiltonian can be written as 

H = H 1 + H 2 -Jd 2 x £ et a eY -/^EE X f-> : < r: y ( (3.2) 

where i, j = 1,2, a, b runs from 1 to M + iV depending on the values taken by the lower 
indices i, j and finally a, (3 run from 1 to n. We consider in the following the replica symmetry 
scheme. Moreover, a recent work on the random bond Potts model has been performed in 



order to test the replica symmetry breaking in [23]. The results so obtained are in favor of 



a non replica symmetry breaking senario. The hamiltonian ( |3.2|) has therefore six coupling 
constants g 1 = g u , 92 = 922, 9i2 = ffai an d &i = A^' 13 , A 2 = A 22 a]b ' fS , A 12 = A 
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^a,a,b,p ^ j g therefore a generalisation of the study of N Ising or Potts models with disorder 



m 0- 



3.2 Beta functions 

As the model we consider has six coupling constants, we have six beta functions. We have 
shown in appendix D how can we obtain them for M = N = 1. Th general case goes along 
the same line but is much more lenghty. Taking directly the replica limit n — > 0, we obtain 



Pgi = Jltl = e i9x + (M-2)g 2 1 -2g l A l + Ng 2 12 

+4A^i - (2M - 5) A l9 \ — (M — 2)gf - Ng l9 j 2 - 2N 9l9l2 A 12 (3.3) 



(3 Al = = e l A 1 - 2A\ + 2(M - l)A l9l + 2N 9l2 A 12 



d Inl 



+2A\ — (M — l)^Aa(2Aa + 9l ) - iV^ 12 A 1 (2A 12 + g l2 ) (3.4) 



Pa2 ^Jltl = ^92 + (N-2)g 2 2 -2g 2 A 2 + Mgl 2 



+4Ajg 2 - (2N - 5)A 2 g 2 2 - (N - 2)g 3 2 - Mg 2 g\ 2 - 2Mg 2 g 12 A 12 (3.5) 



/3a 2 = —4 = £2A 2 - 2A* + 2(N - l)A 2 g 2 + 2Mg 12 A 12 



d Inl 



+2AI — (N — l)g 2 A 2 (2A 2 + g 2 ) - Mg 12 A 2 (2A 12 + g l2 ) (3.6) 
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= = € ~ 1 ^ 1 9i2 + (M-l)g 1 g 12 + (N-l)g 2 g 12 -g 12 (A 1 + A 2 ) 

- l -(M + N- 2)g\ 2 - i(M - l)g l2 gl - l -g l2 g 2 2 + ±(A* + A 2 2 )g 12 (3.7) 
+3A\ 2 g 12 -(M + N- 3)A 12 g 2 12 - (M - l) gi A l9l2 - (N - l)g 2 D 2 g 12 

p Ai2 = ^ll = !l±^A 12 -A 12 (A 11 + A 22 ) + (M-l)^ 1 A 12 + (iV-l)^A 12 



d Inl 



+Mg 12 A 1 + Ng 12 A 2 + A 3 12 + ^A 12 {A 2 ± + A 2 2 ) - -{M - 1)^ 1 A 12 (2A 1 + g x ) 
-~(JV - l))g 2 A 12 (2A 2 + g 2 )- l -{M + N)A 12 g l2 {2A l2 + g 12 ) (3.8) 



Notice that we recover the results of Pujol |L2| by taking in these equations the limit g\ 2 = 
A12 — > 0. We now want to exploit these equations in order to analyse the fixed point 
structure in some peculiar cases and in the general one. 

3.3 Study of the fixed point structure 

We first consider some peculiar values of M and N. The most simple non trivial case consists 
in coupling one gi-state to a q 2 state Potts model under disorder [Tj|. 

3.3.1 Case M = N = 1 

In this case, the beta functions ( |3~3] - |3.8D reduce to 

13912 = dB = e - 1 ^9i2-gi2(A 1 + A 2 ) + ^g l2 (A 2 l + A 2 2 ) + ?>A 2 l2 g 12 + A 12 g 2 l2 

£ A = ^1 = ei A l - 2A\ + 2# 12 A 12 + 2A\ - 2A 1 A 12 g 12 - A x g{ 2 
dim 

/5 Aa = ^f = e 2 A 2 - 2A\ + 2g 12 A X2 + 2A\ - 2A 2 A 12 £ 12 - A 2 g\ 2 (3.9) 
dim 

/3 Al2 = ^ = !i±i^A 12 + ( ?12 (A 1 + A 2 )- A 12 (A X + A 2 ) 
dim 2 



+A 
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~(A? + A 2 ) + A 2 12 - 2A 12 g 12 - g\ 2 



When we consider one Ising model coupled to a 3 or 4-state Potts model, we take the limit 
ei — > and e 2 = e, we find the following fixed points 

g* 12 = ; A* = A 2 = A* 12 = (3.10) 
g\ 2 = 0; A* = 0; A* 2 = \ + \ + o(e 2 ); A* 12 = + o(e) (3.11) 



g\ 2 = 0; A* = 0; A* 2 = + t- + o(e 2 ); A\ 2 = ^ + o(e) (3.12) 
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The first one is trivial and unstable. The second one corresponds to a perfect decoupling 
of the disordered Ising and Potts models (|3.11|) and we obtain a new fixed point mixing a 
priori both models. Let us notice that A* 2 is undetermined at one and two loop calculations 
just enable to compute the first order in e. When one considers a 3— state coupled to a 4— 
state Potts model, we only find the fixed points ( 3.10|) , (|3.11| ). 



In order to study the stability of each of these fixed points, we re-express the systems 
( P-9| ) around the above solutions. Therefore, we write g\2 = g{ 2 + ^9vz] \ = A* + SAi 
and keep only the smallest order in e. We thus obtain a linear system 6X = AX with 
X = (8A±, 8A2, SA12, Sgiz)- All the information concerning the stability is contained in the 
matrix A. We can thus see that ( |3.10| ) and ( |3.12j ) are unstable, and that ( |3.11| ) is stable. 
We have represented in Figure 3, the projection on the flow in the (g±2, A 12 ) plane for two 
different initial conditions (points Aq and Bo). We clearly see that the flow first try to go 
away and then is driven by disorder at the origin corresponding to a perfect decoupling of 
the models. As it has been already noticed in [pl| , such a flow is unusual because it violates 



the c— theorem [£3[] . This example constitutes a new example of first order transition driven 
by randomness in a second order one (see Figure 3) and moreover, the models factorize. 

Another interesting case concerns the limit €\ = t2 = e which corresponds to two Potts 
models. We have therefore g i2 — g, A% — A 12 = A 2 = A and the flow reduces to 

eg - 2gA + g 2 A + AgA' 



d9 . r„-0„\-r-^ i'^' 

d\ogl 



" A -- eA - 2A 2 + 2Ag + 2A 3 - g 2 A - 2A 2 g (3.13) 



dlogl 

We recover the equations introduced by Pujol []T2J. There are three fixed points 

3 = ; A = (3.14) 

the trivial one, and 

g = ; A = i + ^ + 0(e 3 ) (3.15) 
the fixed point corresponding to a decoupling of the models and 



2 2 

g = L + (e 3 ) ; A = ^T^ + 0(e 3 ) (3-16) 



a new non-trivial one mixing both models. The study of the stability of the fixed point 
( |3.16| ) is non-trivial. For the initial conditions g(0) < 0, A(0) > 0, the flow is always driven 
in a strong coupling regime. This is in contradiction with the Imry-Wortis argument |J. For 
g(0) > 0, A(0) > 0, a naive numerical study shows apparently that the flow is also driven 
in a strong coupling regime. Nevertheless, we must not forget that the beta functions and 
therefore the fixed points have been computed order by order in e and the numerical and 



18 



analytical fixed points can not correspond at all in a few cases. This is exactly what we see 
in this case. Therefore, we have to take care with a direct interpretation of the numerical 
flow. Moreover, an analytical study of the fixed point ( |3.16| ) does not enable to conclude 
directly since the coefficient An of the matrix of stability A is zero up to e 2 . A solution 
could be to include higher order terms and to determine the fixed points up to e 3 . Obviously, 
the calculations are probably unfeasable. Nevertheless by assuming the general form of the 
fixed point (|3.16|) up to e 3 , we can show that an does not depend on this form and equals at 



3 

lowest order in e y. By injecting this value in the matrix of stability, we see that the fixed 
point fl3.16|) is in fact stable for g(0) > 0. It proves in fact that fourth order terms (three 
loops one) have at least to be included in a numerical analysis of the flow. 

3.3.2 case N > 1 

We consider in this section the case of M Ising models coupled to N g-state Potts models 
with disorder. We are looking for the fixed points of the system (|3.3|43.8| ). 

First of all, when M = 1 and A^ = A 2 = A^ 2 = 0, we recover the case studied in the 
section 2. The value M = 2 is peculiar since it corresponds to the Ashkin- Teller model. 
Secondly, by considering only two Ising models, the beta functions associated to gi, Ai are 
those of ref. P6| . In this case, the flow can be exactly integrated. When adding the other 
coupling constants, the situation gets quite more complicated. We have done directly a 
simulation of the flow in this case. We also find that disorder makes the models decouple as 
in the previous subsection. 

Otherwise, for general M^2, the fixed points are, up to second order in e, the following 
(the trivial fixed point is still present) 

9*i = , A* = ; 

-e e 2 

gl = + - , A; = if iV > 2; (3.17) 

y2 (N-2) (iV-2) 2 ' 2 ' K ' 

g* 12 = 0, A* 12 = 
g\ = , At = ; 

g * 2 =0 , A* 2 = e - + e -; (3.18) 
^ 2 = 0, A^ 2 = 0or^e 

9*i = , A* = ; 

* e 2 e (3iV-2)e 2 

g 2 = — ,A 2 =~+ ^ , (3.19) 
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= o , A 



12 



or e 

4 



The three different fixed points we have found correspond in fact to the three non-trivial 



fixed points of the "N colored Potts model with disorder [12]. As in the previous subsection, 
we recover for the fixed points ( p. IS 3. IS ) a non trivial value of the disorder mixing a priori 
both models, namely A* 2 = ^e. Nevertheless, it is straightforward to show that this value 
is unstable. We can also prove that the fixed point (|3.18|) is unstable. Therefore, only the 
fixed points ( |3.17| ) and ( p,19| ) (with A^ 2 = 0) are stable. The flow reaches them depending 
on initial conditions (see the discussion in [0). If N ^ 2 and #2(0) < 0, the flow is always 
driven in a strong coupling regime as it was already mentioned in the previous subsection 
(in fact the fixed point ( |3.17| ) goes to infinity). 

Consequently, the infrared physics of M Ising models coupled to iV q— state Potts models 
under weak disorder corresponds to M disordered Ising models plus N disordered Potts 
models. These results can be easily extended to M 3— state Potts models coupled to N 
4— state Potts models. We also find that disorder makes the models decouple. 



4 Conclusions and Discussions 

In this article, we have studied the critical behavior of different coupled minimal models 
under weak disorder. But already in the pure case, we obtain some interesting features. 

By coupling one Ising model to a several Potts models, we get a non-trivial flow of 
renormalisation in a three parameter space (see Figure 1). We have found new tricritical 
points that should correspond to a priori new conformal theories. Moreover, these results 
prove that an energy-energy coupling term is a relevant perturbation which can influence the 
critical behavior. Such a situation might also be encountered in the problem of the XY-Ising 
model. Indeed, a tricritical point close to the fixed point corresponding to a decoupling of 
both models can strongly alter numerical results based notably on finite size scaling. It could 
explain why continuoulsly varying critical exponents are found numerically ]T9 |. Moreover, 



the study of the pure model is also interesting in order to better understand the behavior of 
several coupled conformal field theories. This approach is complementary to the technics of 
exact integrability developped recently in ]17|, |I~6|j . 

We have also seen that the perturbation term can also drive the system in a strong 
coupling regime (it depends on initial conditions) indicating probably a first order transition 
and mass generation. For example, when one considers two different g-state Potts models, 
the flow is always driven in a strong coupling regime. By adding weak randomness, we have 
found that the first order transitions are often transformed in second order one (see Figure 
3) as it should be according to Imry-Wortis argument. Nevertheless, when we consider 
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two coupled Potts models, this argument is violated for some initial conditions in the flow 
of renormalisation. Moreover, by considering the general case of M qi-state Potts models 
coupled to N estate Potts models under weak randomness, we have show that disorder 
often makes the models decouple (exept the above case violating Imry-Wortis argument). 

Recent numerical results of the disordered q > 4-state Potts model have shown a puzzling 
property of the central charge of this model namely c(q — qi x q 2 ) = c(qi) + c(g 2 ) jl5], 0. 
Therefore, we could have expected some relations between the critical behavior of for example 
the disordered 6— state Potts models and one Ising model coupled to 3-state Potts model 
with disorder since their central charge are similar. Nevertheless, the magnetic critical 
exponent (3/v ~ 0.142 found in 0, for the disordered 6— state Potts model appears clearly 
different from the one of the Ising or disorderd 3— state Potts model. Hence, it seems that no 
relations emerge at a perturbative level between two coupled qi and (^-state Potts model and 
a disordered q± x g 2 -state Potts model. Nevertheless, for a large coupling between both Potts 
models, we could expect some relations between the two situations, but it is outside of our 
perturbative analysis. That is why it would be very interesting to investigate numerically 
the critical behavior of several q— state Potts models with disorder especially in the strong 
coupling regime to analyse if we could observe a cross-over phenomena between the two 
regimes. 
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A Beta functions in the pure case 



In this appendix, we compute the beta function of M Ising models (indiced by 1) coupled 
to iV Potts models (indiced by 2). We have three coupling constant g%, g 2 , g% 2 . Some terms 
have already been computed in ||, nevertheless, we recall some details for completeness. 



A.l First order 

A first order correction to g\ comes from the contraction of two {e\E\) terms (plus some 
combinatorial factor corresponding to all of the possible contractions) : 



0^2 ,. M 

J2 e1{x)e\{x)Y,e c x{y)4{y) <?*d 2 y 



\x-y\<r a ^ b=1 c + d 



2 

When b = c ^ a,d this contractions gives : 

A x( r > e i)(^i) 2 / E e^(x)ef(x)d 2 x 

Here, A\(r, ei) (the lower indice corresponds to the first contribution while the upper indice 
indicates first order) is the result of the integral produced by the contraction of two e\ 
operators : 

A\(r,e l )=2(M-2) J < e 1 (x)e 1 (y) > d 2 y = 4tt(M - 2) J (y^J 



\y-x\<r V< r 

-3e a 

3ei 



4vr(M-2)(- 1 (A.l) 



Another first order correction to g\ comes from the contraction of two (£2^2) terms: 

,0 \2 . M AI+N M M+N 



2 

When b = d this contractions gives : 



/ E E E E sl(y)s d 2 (y)d 2 xd 2 y 



\x-y\<r a=1 b=M + 1 c=ld=M+l 



2 X 



with 

dy 



A*(r,e 2 ) = 2N J < e 2 (x)e 2 (y) > d 2 y = 4ttN J ( 



\y-x\<r V<r 
' y— 3e2 ' 



V 



-^Ni— (A.2) 
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We obtain similar correction to #2 by changing (1, M) <-> (2, N). 

There are two contributions to gu coming from the contraction of two (ei£i) or two (£2^2) • 
The integral involved are therefore the same as above and we find 

(Ai(r,e 1 )^ 1 + ^(r ) e 2 )^ 2 °) / ^e^e^x 

J a,c 

with 

A 1 3 (r,e 1 ) = 2(M-l)f^) (A.3) 

Al(r,e 2 ) = 2(iV-l)(^) (A.4) 

A. 2 Second order 

The second order corrections to g^ comes from the contaction of three (ee) terms. We first 
detail the second order terms for g\. 

Contribution to g\ 

• The first contribution comes from the contraction of 2 x 2 E\ operators: 



3! 

\x-y\,\x-z\<r 



r M 

J E ei(^(x)E^(?/)*)E £ i(^(W rf2 ^V^ (a.5) 



The first one corresponds to the case when b = c; d = e; b, d 7^ a, / . We find: 

Al(r,e 1 )(g° 1 f f^e^e^x 

where A\(r, ei) corresponds to the following : 

A?(r,ei)=4(M-2)(M-3) J < e 1 (x)e 1 (y) > < e 1 {y)e 1 (z) > rf 2 ^ (A.6) 

We replace < e{x)\e{y)i > by |x — y\~ 2 ~ 3tl . This integral has been performed in |J and we 
have 

A 2 (r,e 2 ) = 8vr(M-2)(M-3) J j \z\~^\z - l\~^d*z 

y<r ^ ' 

= l 67 r 2 (M-2)(M-3)(^) (A.7) 
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The second term is produced when a = c = e ^ f;b = d ^ f . This will be denoted by : 

Al(r,e 1 )(g° 1 f [^e^e^x 



with: 



A 2 {r, ei) = 4(M - 2) j < e l {x)e l {y)e 1 {z)s x {oo) > < e x {y)e 1 {z) > d 2 yd 2 z (A.8) 

1 3/— x|,|z— x\<r 

The calculation of this integral has been made in 0]. We will nevertheless give the main 
steps because we will have to compute other integrals using similar technics. 
The < £i(x)£i(y)ex(z)Ei(oo) >o term corresponds to the result of £i£\£i — > £\ obtained by 
projecting £\£\£\ over e 1 (oo). We know use the Coulomb gas formulation 

Al{r,e x ) = 4(M-2)Afx 
J < V-^(x)V l2 (y)V 12 (z)V + (u)V l2 (oo) > \y - z\~ 2 -^d 2 yd 2 zd 2 u, (A.9) 

\y-x\,\z-x\<r 

where Af is a normalisation constant which is fixed by taking the limit R — > in the four 
points function. Indeed, we thus have 

< e 1 {0)e 1 (R)e(x)e{y) >= ^ j < V^)V 12 {R)V 12 {x)V + {u)V l2 {y) > d 2 u . 

By taking the limit R — > in both sides, we therefore obtain 



In order to compute ( A.8|) , we first redefine y^y — x; z^z — x. Then, the idea is to pass 



to the "gauge" Z\ = 0, z 2 = z, z 3 — 1, z± — > oo. For arbitrary value of Zi\ i G [1..4], we can 
show using different changes of variables that 

d 2 u < V ai {z 1 )V a2 {z 2 )V a3 {z^)V ai {z i )V + {u) > 



I |4oi02|~ |4aiQ3| |4a 2 c*3 
OC \Zl2\ \Zl3\ \ z 23\ 

|_ |2+4(ai+o 4 )a + +4oi04 1 |4a 2 "++4«2a4 1 ~ |2+4(o3+Q 4 )a + +4a3a4 
\ z 14\ l z 24| l z 34| 

d 2 u\u\ 4aia+ \u - r]\ 4a2a +\u - l| 4a 3«+ ; (A.ll) 

with a\ + a 2 + a 3 + a 4 + a + = 2a as it should be, and n = Z12ZM . Using the above formula 
( |A.11| ) and two trivial changes of variables, we find 



Al(r,ex) = 8n(M-2)Afx (A.12) 

/ / kr 4Ai2 k - ir 4Ai2+4a? Hwi 4a+a ^k-ii 4Q+Qi2 iw-^i 4a+Qi2 rf 2 ^ 



j/<r 
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The second integral is computed in || and recalled in appendix E (eq: ( |E.12| ). It equals 

j— r 4 (— -) 

v 3vr^— 4y. In the z integration, we introduce a new singularity at infinity which need to be 
substracted (see B for more details). We finally obtain 





r 2 1 


V 3e a J 


L 3eiJ 



^(r, ei ) = 87r 2 (M-2) 



•The second type of terms comes from two E\ and two e 2 contractions, namely 



(A.13) 



J Y, e^elix^^e^y^e^e^x^z (A.14) 

11 1 a7^6=l c,d e, f 

las— lf|,|x— z|<r ' ' ' J 

The first contribution comes from the case b = c 7^ a, e; d = f. We denote this term by 

^|(r,6 1 ,c 2 )^^) 3 /5: e ;(x) e ;(x)d a a; 



with 



Al(r, e ls e 2 ) = 8N{M - 2) | < e 1 (x)e 1 (y) > < e 2 (y)e 2 (z) > d 2 zd 2 y (A.15) 

[y-x[,[z-a;|<r 

Using standard manipulations, we find 

„— 3(ei+e2) 



^(r,ei,e 2 ) = 32tt 2 N(M - 2)- 



9eie 2 

The second contribution comes from the case a — c — e;d — f and is denoted by 

;2 



(A.16) 



Al{r,€ 1 ,€ 2 )g 1 ig i3 ) 2 J J2^K(x)d 2 x 



with 



A 2 4 (r, e 1; e 2 ) = AN J < e 1 {x)e 1 {y)e 1 {z)e x {oo) > < e 2 (y)e 2 (z) > (A.17) 

x|,|z— x\<r 

Using the same transformations explained in the calculation of A 2 , we find 



Al{r,e 1 ,e 2 ) = 8ttN J\f x 



(A.18) 



J ( J |^|- 4Ai2 |z-1|- 4A ' 12 +4 ^| m | 4q +^| m -1| 4q + Q12 

y<r ^ ' 



u — z 



4a+ai2 d 2 zd 2 u 



with 4A' 12 = —2 — 3e 2 while 4A 12 = —2 — 3ex. The first integral equals 
second integral is computed in the appendix E (eq: ( |E.16|) ) and equals 



r -3(e 1 +e 2 ) 

3(ei+e 2 ) 



. The 



7T 



V3 



r 4 (-i; 
r2( _ i: 



1 + 



1 1 

3ei 3e 2 



(A.19) 
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As for the calculation of A 2 ,, the z integration in ( |A.18| ) is performed over the whole complex 
plane, while the domain of integration is in fact restricted to the disk \z\ < rjr. This 
introduces a new singularity at infinity. This is equivalent for this integral to add the 
singularity at the origin (see ||). Some easy manipulations show that this extra contribution 
from z — > introduces an integral similar as in ( |A.15| ). We finally obtain 



Al(r,e 1 ,e 2 ) = 16tt 5 



I67H 



_— 3(ei+e2) 

3( ei + e 2 ) 
3( ei + e 2 ) 



1 + 
1 + 



1 1 

3ei 3e 2 
2 

3e7 



+ 16tH 



r — 3(ei+e2) 

9eie 2 



(A.20) 



• There is finally a last contribution to g\ corresponding to the contraction of 2 x 2 e 2 
operators. It comes from 



9l 



M+N 



\ x _ y \ t \ x _ z \ <r a^b=M+l c,d ej 



The contribution comes from the case a = d;b = f. We denote this case by 

The computation of A\{r, ex, e 2 ) is straighforward and we find 

Al(r,e 1 ,e 2 ) = lQir 2 N(N-l) 



9e 



(A.21) 



When we gather all these contributions to g±, we find the equations Q2.9Q . Obviously, by 
exchanging (1,M) <-> (2,N), we obtain the whole contribution to g 2 (see eq. Q2.10|) ). 

Contribution to g 12 



The first contribution to g\ 2 happens with contractions of e\ operators. It comes from 

u 



\x— y\,\x— z\<r 



a^b=l 



We then have two possible contractions. We write the result in the form 



We then follow calculations similar to A\ and A\ (except the symmetry factors) and find 



Al(r,e 1 ,e 2 )=47i 2 (M-r 



-6ei 



36! 



1 + 



2(M-1) 
3ei 



(A.22) 
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We have also the symmetric contribution coming from £ 2 contractions, namely 



A 2 7 (r,e 1 ,e 2 ) = 4n 2 (N-l) 



-6e2 



3e 2 



1 + 



2(N- 1) 
3e 2 



(A.23) 



We then have the contribution coming from the contraction of two e\ and two £ 2 

M M+N 



9i929i2 I E z a i(x)4(x) E el{y)e d 2 {y)Y,e e 1 {z)e{{z)d 2 x<fyd 2 

l„_„l \t-x\^ »#6=1 c/<f=M+l ej 



Z 



\x-y\,\x-z\<r 

By taking a = e 7^ b; c = f 7^ d, we obtain 

a,b 

with 

^(r, ei ,e 2 ) = 167r 2 (iV-l)(M-l 



2 



r -3(ei+£2) 

9eie 2 



(A.24) 



Finally, we have terms coming from (<?i 2 ) (written as (#i 2 ) Ag(r, ei, e 2 ) 

a,b 



(9 



\3 
12 J 



E (*) E (y)^(y) E ^ (^dV 2 - 



|x— 2/|,|x-2|<r a '' 



c,ti 



2,/ 



We have three different possibilities: first contracting two £1 and two £ 2 , second contracting 
three E\ (by the operator algebra E\E\E\ — > £1) and two £ 2 and third contracting three £ 2 
and two £1. The first one leads to a similar contribution as A\. The second and third case 
involve similar integrals as for A\. We find 



A 2 (r, ei ,e 2 ) = 16vr 2 (iV-l)(M-l) 



r -3(ei+e 2 ) 

9eie 2 



+8tt z 



r -3(ei+e2) 

3( ei + e 2 ) 



2(M-1) 2(iV-r 



3ei 



3e 2 



(A.25) 



By gathering all these contributions, we recover the equation (|2.11|) . 



B Renormalisation of rrii in the pure case 

This appendix is devoted to to the computation of the renormalisation of the coupling 
constants associated to the energy operator £j. As in the appendix A, only contraction 
between different £j are involved. 



27 



B.l First order 

There are two contributions at this order to mi. The first one comes from the term 



M M 



9i< J E e* 1 (x)e»(x)'£et(y)<Px<Py 

\x-y\<r a ^ b=1 c=1 

When b = c, we can contract two S\ operators : 

J a 

with 

B\(r, eO = 2(M - 1) j < > d 2 y = -4tt(M - 1) ^^-J (B.l) 

|y-rr|<r 

The second contribution comes from the contraction of an (g\ 2 )eiE 2 with 171262 

r M+N 

9i2m° 2 / E^)*) E eSM^V 

1 1^ a, ft C =M+1 

When 6 = c, we can contract two e 2 operators : 



-3ei 



g° 12 m° 2 B 1 2 (r,e 1 ,e 2 ) f^s^d'y 

J n 



with 

B\(r, ei, e 2 ) = 2N J < e 2 (x)e 2 (y) > d 2 y = -47riV ^ (B.2) 

|j/-x|<r 

This is precisely this term which mix both masses. Obviously, we have by symmetry similar 
contibutions for m 2 . 

B.2 Second order 

We now compute the contributions at second order to m\. The first involves only e 1 operators 

(0)2 , M 

\x-y\,\x-z\<r a ^ b=1 c + d e 

This expression is very similar to the one involved in the calculations of A\ and A\. Using 
the same reasoning, we find 
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with 



S?(r,ei) = 16vr 2 (M- l)(M-2) 



-6ei' 



9e 2 



+ 4?r 2 (M- 1) 



-6ei' 



36! 



3ei 



4tt 2 (M- 1) 



-6ei' 



36! 



1 + 



4M- 6 



36! 



The second contribution comes from 



In this case, we then follow the calculations of A| and A 2 . We find 



(^ 2 ) 2 m?5 2 2 (r,e!,6 2 ) / 5>?(s)rf 



x 



with 



5 2 (r,ei,e 2 ) = l6n 2 N(M - 1) 



' — 3(ei+£2) ~ 



96i6 2 



8vr 2 iV 



r 



-3(ei+e 2 ) 



,3(ei + e 2 ) 



1 + 



36! 



4n 2 N 



r -3(ei+£2) 



1 + 



2M 2(M-1) 



3ei 



3e 2 



,3( ei + e 2 ) 

We have finally the contributions coming from m" at second order, namely 

M AI+N 



| as— |a— z|<r 



and 



M+7V 



M+7V 



E 4W4WE^(2/)*) E e e 2 (z)d 2 xd 2 yd 2 ; 

\x-y\,\x-z\<r <^b=M+l c,d e+M+1 



Taking b = c; d = e in the first one and a = e; b = d in the second one we obtain 
(^X^3 2 ('",6i,e 2 )+^ 2 X5l(r,e 1 ,e 2 )) f 5>?(x)d 2 x 



(B.3) 



(B.4) 



with 



B 2 {r,e 1 ,e 2 ) = 8n 2 N(M-l) 



Bl(r,e 1 ,e 2 ) = 8n 2 N(N - 1) 



3(ei+£2) ~ 

v 9eie 2 



9e 2 



(B.5) 
(B.6) 



By collecting all these contributions, we obtain the equations Q2.16p and by symmetry the 
equation ( |2.17| ). 
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C Renormalisation of hi in the pure case 



We now compute the renormalisation of the coupling constants hi(r) associated to the spin 
operators cr.j (<7i the Ising spin operator, <ji a Potts spin operator). We restrict ourself to 
the most interesting case, namely M — 1 Ising model coupled to iV Potts models. Here, 
the renormalisation of the spin operators is induced by the interaction term (see eq: ()) 
containing only energy operators. Therefore, we have a mixture of <7j with £j. The first order 
term gives no contribution because the OPE of two e operators does not involve a operators. 
The first non-zero contribution is at second order. 



C.l Second order 

Here, we compute the renormalisation of a Potts spin operator with the interaction term 



YjYj9ij I d 2 z £ ii z ) £ j{ z )- There are two non-zero contributions. The first one comes from 



c,d i,j 



(0?) 



0\2 







J £ e a x {x)e\{x) £ e?(y)ej(y) £ a\{z)d 2 xd 2 yd 2 z {C.l) 

\v.— 7.\<rr a + h c + d e 



\x—y\,\x— z\<r 



A similar integral has been encountered in f5j. When a = c = e^b = d the product 
(j\e\E\ contains a a operator. As usual, we project the product <J\e\E\ on a 01(00) operator, 
obtaining therefore 



„ M 

hMfC?{r,e u e 2 ) / E^iW^ 2 - 

J a=l 



where 

C{ 2) (r,ei,e 2 ) = 2(iV - 1) J < a 1 {x)e l {y)e 1 {z)a 1 {oo) > < e x (y)e x (z) > d 2 yd 2 z 

\y—x\,\z—x\<r 

(C.2) 

Using the Coulomb gas formulation, this integral reads 

4(M - 1) A/V / (-^)x (C.3) 

y<r ' 

J \ z f a ^ a —\ z - l|-4Ai2+4a2 2 | u |4 a+ a^ 3T | u _ X |4a + «i2 | ^ _ ^| 4a+a 12 

where the normalisation Af is the same as the one computed in the appendix A. The second 
non-trivial has been computed in H (see also the equation (E.23|) for such integrals) and we 
find 



C 1 (2) (r,e 1 ,e 2 ) = (M-l)7r 2 



-6ei 



l + 4 f2 A ) g=lg(l) 

1 + 3 (2 - Al) r 2 H)P(-i; 



(C.4) 
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The other one comes from 

i? / £ el{x)e b 2 {x) £ el(y)e d 2 (y) £ of (^d 2 ^^ 



(C.5) 



| as— z|<r 



0,6 



c,<2 



We take a = c = e, b = d and project on the 01(00) operator and thus obtain 



„ M 

%i° 2 ) 2 ^ 2) (r,£i,e 2 ) /£^)rf 2 - 



a=l 

with 

C| 2) (r,ei,e 2 ) = 2iV / < <n(a;)ei(i/)ei(«)«Ti(oo) > < ^(y)^*) >o dV 2 ^ (C.6) 

|j/-a:|,|*-a![<r 

Using the Coulomb gas formalism, we obtain an integral very similar to ( |C.3|) 

dy \ 



4N Afn 



y<r 



l+3ei+3e 2 



(C.7) 



/ |^| 4ai2Q ^|2 - l|- 4A i2+ 4 °? 2 | u | 4a + a ^T| u - \\ Aa +^\ U - z\ Aa + aw d 2 zd 2 U 



with 4A' 12 = —2 — 3e2 whereas 4Ai2 = —2 — 3ei. This integral has been computed in the 
appendix E (see ( |E.27| )) and yields 



{2) 3N 2 r-^j 

69 (r, ei,e 2 = vr — 

2 V ' ' 7 2 3 ei + e 2 



e 2 + ei(l + -{2 - AiJ r2 ,_i. r2 ,_i, 



Obviously, we get similar contributions for h® when we substitute (1,M) <-> (2,N). 



C.2 third order 

At third order, we have to do the contraction of threee (ee) operators with one a operator. 
We obtain three main contributions. 

•The first one does not mix the E\ and e 2 operators and has been computed in |J. It 
reads as 



Ml 3 

3! 



h\ J £ el{x)e\{x) £ ef (y)e?(y) £ ef(x)e{(x) £ a^xd^zd 

\ x -y\,\ x -z\,\x-u\<r a ^ b 9 



U 



Here contributions come from the following contractions a — c — g;b — e;d — f and 
produce : 



(g° i rh° 1 C!(r,e 1 ) £<(^ 2 " 



o=l 
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where 

C?(r,ei) = 4(M-l)(M-2)x (C.9) 
y < cri(x)£i(2/)£i(2;)cri(oo) > < £iG/)ei(w) >o< e 1 (w)e 1 (z) >o d 2 yd 2 zd 2 w 

\y—x\,\z—x\,\w—x\<r 

and leads in the Coulomb gas representation to 

( dy \ 



(M - 1)(M - 2)Afvr / \w\~ 2 ~ 3tl \w - l\~ 2 ~ 3ei d 2 



w 



y<r 



In order to factorize the y and u> integrations, we have used the formula ([A. 11 ) and three 



trivial changes of variables. The first two integrations will then produce ^|jgr~^ while the 
results of the z and u integrations are given in H and in appendix E ( |E.30| ). Therefore, the 
final result is 



C{ (r, ei) = -12(M - 1)(M - 2)vr 
» The second contribution reads as 



3 



' -9e-l' 



9e- 1 



r 2 (-I)r 2 (|) 



1 + 9 (2 Al) r2(-i)r 2 (-i)j 



(C.ll) 



Taking for example a = c = g; b = e; d = f , we obtain 

. M 

J 0=1 

where 

C 2 3 (r,ei) = 8N(M-l)x (C.12) 
/" < o"i(x)£i(2/)£i(2;)o-i(oo) > < ei(2/)ei(w) > < £2(^)^2(2;) >o d 2 yd 2 zd 2 w 

\y—x\,\z—x\,\w—x\<r 

and leads in the Coulomb gas representation to 



16JV(M - 1)A/V / |w 



| \w\- 2 ' 3 ^\w-l\- 2 -^d 2 w J ( 



^l+6ei+3e 2 
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The third integral is computed in appendix E (see ( |E.33| )) and the final result is 



C|(r,ei,e 2 ) = -8N(M - (g° 12 ) 2 g° 



-6ei— 3e 2 r 



6e 1 + 3e 2 



1,, , e 2 . 



(l + ^)(l + -(2-A 1 )^) + -(l + -) 
e 2 o z ex J 

(C.14) 



The third contribution reads as 





12/ ^OlO 



3! 



9lK 



M+N 



E e a x {x)e\{x) E e?(y)ef (y) E ^(x)ej(x) E af(z)d 2 xd 2 yd 2 zd 2 u 



c.d 



\x-y\,\x-z\,\x-u\<r a ' b c ' d e^f=M+l 

Taking for example a — c — g;b — f;d — e, we obtain 

M 



M 2 9 2C$(r,e lt e 2 ) / E 



2 



0=1 



where 



C|( r , ei ) = 8iV(iV-l)x (C.15) 
< oi{x)si(y)ei(z)ai(oo) > < E2(y)e 2 (w) > < ^(w)^^) >o d 2 yd 2 zd 2 w 



\y— x\,\z— x\,\w— x\<r 

and leads in the Coulomb gas representation to 



16iV(iV - / \w 



|-2-3e a 



W — 1 



-2-3e 2 „r2 



ciy 



.l+3ei+6e 2 



1 ,4ai 2 a r 



— I ^ _ l|-8A' 12 +4a? 2 +2| lt |4a + a^ 5=T | /ii _ 1 |4a+a 12 | u _ ^ 1 4o+«i2 ^2^2^ _ ^ ^ 



with 4A' 12 = — 2 — 3e 2 . The third integral is computed in appendix E (see ( K.3flp ) and the 
final result is 



Cl(r, e 1; e 2 ) = -8N(N - 1)^)^) 



-3ei— 6e 2 



3ei + 6e 2 



l + ^(l + ^(2-A 1 )^) 



(C.17) 



When we gather all these second and third contributions, we find the renormalisations 
of h°i fl2.19U2.20Q . 



D Beta function in the disordered case 

In this appendix, we only give the renormalisation associated to the coupling constants of 
one Ising model plus one Potts model with disorder, namely g i2 , A 1; A 2 , A 12 . We do not give 
the details of all contributions but only the final results since the integrals we compute are 
similar to appendix A. The only difference lays in symmetry factors. The results extend to 
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M Ising models coupled to N Potts models with disorder, nevertheless there are much more 
terms that will simplify when we use the renormalised coupling constants. In the following 
equations, n represents the replica number. 



Ai(r) = r 



-3ei 



A o _ 4n{n _ 2 )(A?) 2 — - - 4vm(A? 2 ) 2 — - - 87r<&A° 2 — - - 

061 Ot 2 3E 2 

+ 4(n - 2)(ra - 3)(A°) 3 J n + 4(n - 2) (A?) 3 Jn 

+ 8(n - 2) 2 A?(A° 2 ) 2 J 12 + 4(n - 1)A°(A° 2 ) 2 Ji 2 (D.l) 
+ 16(n - 2)A?A? 2 (A? 2 + ^ 2 )J 12 + 4(A? 2 + <?° 2 ) 2 A?J 12 

+ 4 ( (2n 2 - 3n + 3)(A? 2 ) 2 A° + 2(n - 2)(«& + A? 2 )A° 2 A° + (g° 12 + A° 2 ) 2 A°) J 22 



A 2 (r) 



-3e 2 



A° - 47r(n 



2)(A°) 2 - 



-3e 2 



3e 2 



4vm(A 



-3ei 



12; 



36! 
. 0\3 



-3ei 



3ei 



+ 4(n - 2)(n - 3)(A^)% 2 + 4(n - 2)(A 2 , ) 3 J 22 
+ 8(n - 2) 2 A°(A° 2 ) 2 J 21 + 4(n - 1)A°(A° 2 ) 2 .7 21 (D.2) 
+ 16(n - 2)A°A? 2 (A? 2 + g° l2 )I 21 + 4(A? 2 + g° l2 ) 2 A° 2 J 21 

+ 4 ((2n 2 - 3n + 3)(A? 2 ) 2 A? + 2(n - 2)(<7? 2 + A? 2 )A? 2 A? + (<& + A? 2 ) 2 A?) J n 



(A 12 + # 12 )(r) = 



-f (ei+e 2 ) 



f 2 



A 12 (r) 



3ei j,— 3e2 

(A° 2 + - 4vr(n - 1)A?A? 2 ^— - An(n - 1)A°A° 



36! 



i 12' 



3e? 



+4(n - l)(n - 2)(A?) 2 A° 2 J n + 4(n - l)(n - 2)(A°) 2 A° 2 X 
+2(n - 1)(A?) 2 (A? 2 + <?? 2 ) Jii + 2(n - 1)(A°) 2 (A? 2 + «&) 
+4(n -i)(„- 2)A°A°A° 2 J 12 + 4(n - 1)A?A°(A? 2 + ^° 2 )J 12 
+4(n -i)(„- 2)(A? 2 ) 3 J 12 + 4(n - 1)(A? 2 ) 2 (A? 2 + ^ 2 )J 12 
+2(n - 1)(A? 2 ) 2 (A? 2 + ^ 2 )(Ji 2 + J 21 ) 



22 



(D.3) 



AO 

12 



-4ttA?(A? 2 + ^ 2 ). 



A° 2 - 4n(n - 2)A°A 

„— 3ei 



3ei ™— 3e2 

- 4vr(n - 2)A°A° 



3ei 



x 12' 



3e 5 



36! 



4ttA°(A? 2 + ^ 2 ). 



-3e 2 



3e 2 



+4(n - 2) 2 (A?) 2 A? 2 J n + 4(n - 2)(A?) 2 (A? 2 + ^ 2 )J n + 2(n - l)(A°) 2 A° 2t 7 u 
+4(n - 2) 2 (A°) 2 A? 2 J 22 + 4(n - 2)(A°) 2 (A? 2 + g° 12 )l 22 + 2(n - 1)(A°) 2 A? 2 J" 22 
+4(n - 1) 2 A°A°A° 2 X 12 + 4(n - 2) A°A°<?° 2 X 12 (D.4) 
+4(n - 2)(n - 3)(A? 2 ) 3 J 12 + 2(n - 2)(A° 2 ) 3 (J 12 + J 21 ) 
+12(n - 2)(A? 2 ) 2 (A? 2 + £° 2 )J 12 + 4(A? 2 + ^ 2 ) 2 A? 2 J 12 
+2(A? 2 + ^ 2 ) 2 A? 2 (J 12 + l 7 21 ) 
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In this system of equations, we have defined for readability the following integrals 



Dei 



11 



4tt z 



9e 2 



j ^22 



4^ 



9e2 



11 



2tt 2 



J12 = 4vr' 



36! 

r -3(ei+e2) 

3( ei + e 2 ) 



12 — 


x 21 






2 " 




1 + 






3^_ 






2 " 




1 + 


3^- 





4tt 2 



-3(ei+e 2 ) 

9eie 2 



22 



21 



2vr 2 



4tH 



,r 



-6e2 



3e 2 



3e 2 

r -3(ei+e2) r 



1 + 



3e 2 



(D.5) 

(D.6) 
(D.7) 
(D.8) 



3(ei + e 2 ) 

The calculations of the last four integrals is recalled in the appendix E. The last step consists 
in deriving the beta function, r being the scale variable. 



E Computation of integrals 

In this appendix, we detail the calculation of some integrals we have used in the previous 
appendices. These integrals take the general form 

I = J \ x \ 2a \x - l\ 2b \y\ 2a '\y - l\ 2b '\x - y\^d 2 xd 2 y (E.l) 

where the integration is performed over the whole complex plane. Some general techniques 
to compute this kind of integrals have been shown in ref . p4], [27]. We just recall the main 
steps in this appendix (see also ref. |J for the method). This integral can be decomposed 
into holomorphic and antiholomorphic parts : 

I = s(b)s(b') [j^Jf + 4J2] + s(b)s(b' + 2g)J+J 2 + s(b + 2g)s(b')J+J^ (E.2) 

where s(x) corresponds to sin{irx) and 

J+ = J ( a , 6, a, 6', g) ; J 2 + = J (b, a, b', a, g) 

J- = J(b,-2-a-b-2g,b',-2-a' - b' -2g,g) (E.3) 
J- = J(-2-a-b-2g,b, -2 - a' -b'-2g,b',g) 

Here, we used the notation 
1 1 

J(a,b,a',b',g) = J du J dv u a+a ' +29+1 (l - u) V (1 - vf 9 {\ - uvf (E.4) 


r(2 + a + a' + 2g)T(l + b)T(l + a')T(l + 2g) ~ (-b') k {2 + a + a' + 2g) k {l + a') k 



E 



r(3 + a + a 1 + b + 2g)T(2 + a' + 2g) ^ fc!(3 + a + a' + b + 2g) k (2 + a' + 2g) k 
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and 

(a)k = a(a + l)...(a + k — 1) 

The J integrals appearing in (|E.2| ) are not all independent. Using contour deformation of 
integrals it can be shown that we have the following relations: 

s(2g + a + b)Jr + s(a + b)J 2 = ( s ( a ')J+ + s(2g + a') J 2 + ) (E.5) 

s(2g + a' + b')J 2 + s{a! + b')J{ = s(2 ^ + h) + s(2g + a)J+) (E.6) 

With these formulas, the integrals appearing in the appendices A,C can be calculated as a 
power series of ei and e 2 . 

• Let us detail the integral referred as ( |A. 12|) . This integral has been computed in || 
but we recall it since we find some errors in the intermediary steps of ref . jS| . 
The coefficients corresponding to that case are 

a = -2a' = 2 + 3 ei ; b' = ~ - e x ; b = 2g = -~e 1 (E.7) 
Substituting these values in QE.2Q , we obtain at the second order in e\ 

T - 7(j 2 + -j 1 + )(Jr-^)-^^(^(2jr + j 2 -) + j 2 + (jr + 2j 2 - 

vr 2 e? 



-A (4(J+Jr + J+J 2 ) - 19(J+J 2 + J+Jf)) (E.8) 



Now we write the J i as functions of the Jf with the help of relations (|E.5|JE.6| ) 

J- = _j t + 2 U_^_tf\ J} (E . 9) 



4 



J2 



' r 19 , 2 \ . / 3vr\/3ei 33vr 2 e 2 \ . , 

1 + vrv^ei + y7r 2 e 2 j J+ +2 1 + — ^ + — ^ j J 2 + (E.10) 

The equation (|E.S| ) simplifies to 



/ = 3y/3ne 1 J?J+ - ^7r 2 e 2 (J+) 2 (E.li; 



The Ji" and J 2 can now be computed explicitly using formulas ( |E.3|JE.4| ). We obtain : 



4 r 2 f-i) 

O Q p2/_ 1 \ 

4 = g(2 +TrV ^ ) ^(I) +0{e] 
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The final result is therefore (at lowest order in ei) 



r 2 (-i: 



• We now compute the integral ([A.18Q . The coefficients are now 
a = -2a' = 2 + 3e 1 ; b' = -~ + ^ - 3e 2 ) ; b = 2g = - 
Substituting these values in ( |E.2| ), we obtain at the second order in e 



(E.12) 



(E.13) 



7T61 



f^i (19( J+J{ + J+J 2 -) - 19J+J 2 - - 43J+J 1 - 

(5(j+jr + j+j 2 -) + 5j+j 2 - + j+jf 



Ji + J<z J2 5*^2 "^1 



Jl J\ + J% J% + J\ ^2 ^2^1 



+ 



32 

37r 2 eie 2 
16 - 
27tt^| 

32 



J\ J\ + </ 2 </ 2 J\ J2 ^2^1 



(E.14) 



Now we write the Jj as functions of the J 4 + with the help of relations (|E.5| , [E.6|) and the 
equation (|E.14|) simplifies after tedious algebra to 



I = 3V3n^J+J^ 



3V3 



-ne 



£2 4 e 2 ei 2 e 2 

The Ji" and J2 1 " are now be computed explicitly using formulas ( |E.3|JE.4p . We find : 



(E.15) 



Jt 



4 r 2 (-|; 

9e a r(i) 



1 + (~ + + o(e*) 



2 ( 3 +7r ^3 } £!(i) +0( 

9 V 2 ; r I v 



The calculation of the sums 



is quite heavy. We refer to ref. |27[ where some tricks 



concerning their manipulations are given. 
The final result is therefore 

r 4 r-i' 



V3 



71 



F»(-§) 



1 1 
1 + 

3ei 3e 2 



(E.16) 



Note that we recover the result (E. 12|) when ei = e 2 as it should be. During the calculation, 
we have checked that the terms in ^ simplifiy each other in order to ensure renormalisability. 



37 



6 

Nevertheless, we are obliged to admit the result for terms in ^. On the other hand, singurities 
in — are necessary in order to compensate the ones coming from one loop calculations. 



We now compute the integral ( p.3| ). Here, we have : 



a = h Aiei = — 2a : 

3 

Substituting these values, we find 

i = \(4-Jf){j;-j;)-^ 



b = 2g = b' - 1 



(j+(2jr + j 2 ) + j 2 + (jf + 2j 2 



Then, relating the Jj to the Jf, 

J+ = -7r^3ei(2- Ai)Jf - (l-7r^3ei(2- Ai))J 2 " 
J 2 + = -vr^ ei (2-A 1 )Jr-(l + 7r^ ei (l + A 1 ))J 2 - 



we obtain for J : 



/ = r j 2 ) 2 + ^ (2 ~ Al)7rei (j 2 - - j r ) 2 + o(e?; 



2 v ' J 2 
These J~ can again be computed with the help of ([E.3|JE.4| ) : 



J i = J 2 + 



r(-f)r(i)_ir-(-§) , r(-i)r(j) 



+ 



r(-|) ~ 2 r(-l) ^ r(-| 



+ O(ei) 



Therefore : 



_ 3V37rei T 4 (-| 



r 2 (-4) 



4,- - - r 2 (-|)r 2 (i; 



1 + 3 (2 " Al)7rei P(-|)P(4; 



(E.17) 



(E.18) 

(E.19) 
(E.20) 

(E.21) 



(E.22) 



(E.23) 



We now compute the integral ( |C.7| ). Here, we have now : 



a = h XiSi = —2a' ; b = 2g 

3 

Substituting these values, we find 

3 



ei ; 6' = ~ + ^(6i-3e 2 ) 



- j+)(jr - j 2 ) 

a F^Xttiz^ ^ ^ ^ 

J I J I + J 2 J 2 ~ <A ^ 2 + 5J 2 J i 



(j\ J\ + ^2^2 + ^1^2 ^2^1 



(E.24) 



(E.25) 



38 



The relations between the J i to the Jf reamin the same as above and we obtain for / : 
I = 3V3n{e 1 + e 2 ) / _y | y/3(2 - Ai)tt£i / . . x2 



J 2 J x 



O(el) 



(E.26) 



The values of the J i remain the same as ( |E.22j ) and thus 



_ r 4 (-| 



16 r 2 (- 



r 2 (-f)r 2 ^' 



£2 + ei l 1 + 3 (2-Al) r 2 (-i)P(-|), 



(E.27) 



We now calculate the integral ( U.10Q . We have the following coefficients: 



1 



■ 3 + A £l 



-2a! ; b = 2g 
With similar manipulations, the integral reads 



b' 



1 5 

ei 

3 2 



9v/ ^ ei 2 -) 2 + ^ (2 - Al)7rei (j 2 - - j r ) 2 + 0(e 2 ) 



and the final result is 

T 9V37r ei r 4 (-±) 



16 r 2 (-|) 



8, . r 2 (-|)r 2 (i) 
1 + 9( 2 - A ^ r 2 H)r 2 H) 



+ o(e 2 ; 



(E.28) 



(E.29) 



(E.30) 



We now compute the integral ( |C.13| ). We have: 



a = — - + Aei = —2a' 



With similar manipulations, the integral reads 

T 3 v / 3vr(2e 1 + e 2 ) /_x2 ^3(2 - Ai)7rei 



Jo + 



J 2 --Jf +0(e 2 



and the final result is 
3v37rei T 



4f_ I x 

3' 



16 r 2 (-f) 



/ 4 r 2 ('-^)r 2 (' 1 -) 

e 2 + 2^1 + -(2-A 1 K l r2( ^3^ ) 



+ 0(e 2 



(E.31) 



(E.32) 



(E.33) 



Concerning the integral (|C.16|) , we have the following coefficients: 



-- + M 



-2d 



b' 



1 , ei 



-3e 5 



With similar manipulations, the integral reads 

r 3 v / 3vr(ei + 2e 2 ) / _n2 >/3(2 - Ai)^ 



4 



and the final result is 

_ _ 3^6! r 4 (-|; 



16 r 2 (- 



j. 



/ 8, . r 2 (-|)r 2 (i) 
2£2 + H 1 + 3 (2 - Al) " £l r 2 (-|)r 2 H) 



J--J-) +0(el 



+ 0(e 2 



(E.34) 



(E.35) 



(E.36) 
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CAPTIONS 



FIG.l : The projection of the flow in the (gu, 92) plane for one Ising model coupled to three 
Potts models. Two new tricritical points denoted by Ti, T 2 are found. 



FIG. 2 : The projection of the flow in the (^1,^12) plane for two Ising models coupled to a 
Potts model. The point P separates the fixed point line g 1 in stable and unstable parts. 



FIG. 3 : The projection of the flow in the (<?i2,Ai 2 ) plane for one Ising model coupled to 
one Potts model under weak disorder. We clearly see that disorder prevents the flow from 
running in a strong coupling regime and therefore makes the models decouple. 
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Figure 1 
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Figure 2 




0.05 0.1 0.15 0.2 



§12 



Figure 3 



45 



